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We calculate the nucleon self-energies in nuclear matter in the QCD sum rules 
approach, taking into account the contributions of the four-quark condensates. We 

o : 

■ analyze the dependence of the results on the model employed for the calculation of 

the condensates and demonstrate that the relativistic character of the models is im- 
portant. The condensates are calculated with inclusion of the most important terms 

in 

CN ' beyond the gas approximation. This corresponds to inclusion of the two-body nu- 

cleon forces and of the most important three-body forces. The results are consistent 
with the convergence of the operator product expansion. The density dependence of 

O 

the nucleon self-energies is obtained. The results are consistent with those obtained 

G,; 

by the standard nuclear physics methods, thus inspiring further development of the 
. approach. 

o 

^ ' I. INTRODUCTION 

in 
o _ 

Calculation of the density dependence of the nucleon self energy is one of the most important 
problems of nowadays nuclear physics . The standard nuclear physics calculations carried 
out in various approaches j^, S] require certain phenomenological parameters. The QCD sum 
rules approach to description of the nucleons in nuclear matter [4|] enables to present the nucleon 



self-energies in terms of the in- medium QCD condensates |5j. 

The QCD sum rules approach stimulated investigation of the in-medium four-quark conden- 



sates. It was demonstrated [6] that the influence of the four-quark condensates may appear to 
be large. However the calculations require certain model assumption on the quark structure of 
nucleus. 

In the contribution of the scalar condensate to the QCD sum rules was obtained by using 
the Nambu-Jona-Lasinio model js]. A larger number of condensates was obtained in [s| in 



framework of the Perturbative Chiral Quark Model (PCQM) 10, [ill]- Further studies showed 
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13| that a more detailed analysis is needed. Also, the status of the factorization assumption 



5] should be clarified. 

In the present paper we consider the four-quark condensates in the class of models, where 
the nucleon is viewed as a system of the valence quarks and the pions, which are treated 
perturbatively. The sea quarks are assumed to be contained in the pions. We focus on the role 
of the four-quark condensates in the QCD sum rules analysis of the nucleons in nuclear matter. 

The QCD sum rules approach to investigation of hadrons in vacuum was suggested in [l4| for 



the mesons and was extended for nucleons in 



15| (see also 16] ) . The aim of the approach is to 



express hadron parameters in terms of the QCD condensates. It was extended for calculations 
of the nucleon self-energies in nuclear matter (see 17| and references therein). The method is 
based on dispersion relations for the function IT m (g 2 ) usually referred to as the "polarization 
operator" which describes the time-space propagation of the system with the four-momentum q 
and the quantum numbers of the nucleon in the matter. The function U m (q 2 ) is determined by 
the local quark operator j(x). The choice of j(x) is not unique. We shall use the form suggested 
in [3] (see Appendix A). We shall focus on the contributions of the four-quark condensates to 
the function II m (g 2 ). 

The calculations are carried out in the gas approximation and with inclusion of the most 
important nonlinear terms. This corresponds to inclusion of the two-body forces and of the 
three-body forces in the mean-field approximation. 

The unknowns of the sum rules equations are the nucleon vector self-energy Sy, the effective 
mass m* (or the scalar self-energy m* — m), the nucleon residue and the effective continuum 
threshold W^. The sum rules equations connect these parameters to the QCD condensates at 
finite density. Putting equal to its vacuum value we find explicit approximate equations 
which express the nucleon self-energies in terms of the QCD condensates. 

There are three types of contributions to the four-quark expectation values. In the first one 
only the valence quarks are involved. The second one includes only the sea quarks. There 
are also the interference contributions, where both sea and valence quarks participate. The 
contribution of the valence quarks can be written in a way, which enables to separate special 
factorized configurations in which two quark operators act inside the nucleon while two other 
ones act on the QCD vacuum. 

It was shown in [3| that the contributions involving only the sea quarks cancel. Turning 
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to the valence quarks we find the factorized terms and also the internal terms in which all 
four operators act inside the nucleon. In the interference terms two quark operators act on the 
constituent quark, while the other two act on pion or connect the pion with vacuum by the 
PCAC relations. 

The four-quark condensates have dimension d — 6. The factorized terms are proportional to 
the product (0|gg|0) • (N\qT x q\N) va i ■ p (T x = 1, 7o) of the vacuum expectation value (0|gg|0), 
contribution of the valence quarks to the nucleon expectation value (N\qTxq\N) va i and the 
nucleon density p. The internal terms depend on the characteristic nucleon size R. They scale 
as p/R 3 . Since (0|gg|0) « (-241 MeV) 3 , while R « 0.6 fin = (328 MeV)- 1 , these two values are 
of the same order. 

We demonstrate that in the scalar channel the factorized terms determine about 80 percent 
of the total contribution. The internal contributions are as important as the interference terms. 
In contrary, in the vector channel the nonfactorized terms play the crucial role. There is a large 
cancellation between the factorized and internal terms. 

The results are consistent with the assumption of the convergence of the power series in 1/q 2 
for the polarization operator. 

We show that in the gas approximation inclusion of the four-quark condensates subtract 
more than 100 MeV from the vector self-energy Ey and from the nucleon effective mass m* at 
the saturation value of density. The effect becomes weaker after inclusion of the main three- 
body terms, in which the two pairs of the quark operators act on two different nucleons. The 
quantitative results depend on the value of the nucleon matrix element = (N\uu + dd\N). 
Anyway, at the phenomenological saturation point the matter appeared to be bound. 

In is instructive to compare the results with those, obtained in framework of the nonrela- 
tivistic quark model. The results differ strongly from those found in the PCQM. In particular, 
the matter appeared to be bound only if the value of kn is large enough. 

Our approach enabled to find the density dependence of the nucleon self-energies. The 
reasonable results of the paper stimulate further calculations with more complete inclusion of 
the three-body terms. 
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II. GENERAL EQUATIONS 
A. Expression for the four-quark condensates 

We calculate the expectations values of the color-antisymmetric operators 

T x f = : (q ha T x q ha ' • " <W<W)) : • (1) 

Here the colon signs denote the normal ordering of the quark operators qf , /i j2 stand for the 
quark flavors, while a, a', b, b' are the color indices. The basic 4x4 matrices T x ,y acting on the 
Lorentz indices of the quark operators are 

T x = I, T x = 75 , r x = 7m , = 7^75 , 

Tx = 2^7^-7^) = <V» ( 2 ) 

describing the scalar, pseudoscalar, vector, axial and tensor cases correspondingly. The opera- 
tors defined by Eq. (1) can be written also as 

T S x f = : • q h V x I c q h • g /2 r y / c g /2 - ^ q h V x \ a q h • g*IYA a g*) : , (3) 

where A a are the basic SU(3) matrices, J c is the unit 3x3 color matrix. We shall omit J c 
further. 

B. Sum rules 

L Inclusion of the four- quark condensates 



Considering nuclear matter as a system of 
vector 

P = 

In the rest frame of the matter 

P, = 

with po ?s m. 



A nucleons with momenta pi we introduce the 
A ' 
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The sum rules are based on the dispersion relations 

for the components (q 2 ) of the polarization operator 

IT m (g 2 ) = i [ d*xe^{M\j(x)j(0)\M) = qW m {q 2 ) + PTl P m {q 2 ) + IU'jq 2 ). 



The left-hand side (LHS) of Eq. is presented as a power series of q 2 . This is known as the 
operator product expansion (OPE). The functions Imll^/c 2 ) on the right-hand side (RHS) of 
Eq. (4) are approximated by the "pole+continuum" model in which the lowest pole is written 
exactly, while the higher states are approximated by continuum. 

The Borel transformed QCD sum rules in nuclear matter take the form 

a(M 2 , W 2 n ) = A^exp (-m 2 m /M 2 ); C J (M 2 , W 2 m ) = m'A^exp (-m 2 m /M 2 ); (5) 

C P (M 2 , W 2 m ) = -H v \ 2 m exp (-m 2 m /M 2 ) . 

Here D are the Borel transforms of the LHS of Eq. (4) with the continuum contributions 
transferred to the LHS, m m is the new position of the nucleon pole [17]. These equations are 
expected to determine the nucleon effective mass m* and the vector self-energy Sy, the residue 
at the nucleon pole X 2 m and the effective continuum threshold W^. The explicit form of the 
LHS will be presented later. 

The LHS of Eqs. (5) which includes the in-medium condensates can be written as (we omit 
the dependence on W^) 

C\M 2 ) = Y,MM 2 ); C I (M 2 ) = J2Bn(M 2 ); C p (M 2 ) = £ P n (M 2 ) . (6) 

n n n 

Here the lower indices n denote the dimensions of the condensates. Inclusion of the higher 
dimensions corresponds to inclusion of the higher terms in the 1/g 2 expansion of the LHS 
of Eqs. (4). Each term on the RHS of Eqs. ([6]) is the Borel transform of the corresponding 
contribution to the polarization operator W(q 2 ) multiplied by the factor 32tt a . 

Earlier Eqs. ([5]) were solved for the contributions with n < 6 with inclusion of the radiative 
corrections of the order a s . Now we include also the terms with n — 6, which are mainly the 
four-quark condensates. 
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The contribution of the four-quark condensates to the polarization operator can be written 

as 

(n m ) 4g = (n ) 4q + \j2{t iXY H X Y + T XY R XY ' ) ), (7) 

^ X,Y 

where (n )4 g is the vacuum term, with 

H XY = (M\T X U Y \M) - (0\T X U Y \0); R X y = (M\T x d Y \M) - (0|T^|0) , (8) 
while (see Appendix A) 

V XY = ^ Tr (rT x ^T Y ) 7s7 i W5 , (9) 
lo 



and 



r XY = 6 -f Tr (r^) ^ la Y Y lpl5 . (10) 



Here 9x — +1, if T x has a vector or tensor structure, while 9x = — 1 in other cases. 

Here we omitted the Lorentz indices of the matrices Fx, p XY and t xy . Note the way, 
matrices T x with the lower indices in Eqs. (l)-(3), are related to those with upper indices in 
Eqs. ([9]) and (TlO|) . In the scalar and pseudoscalar cases T x = T x , T x = 7 M and T x = a^ v in 
vector and tensor cases, T x = 757^ in the axial case. 

The traces on the RHS of Eq. (|9]), describing contributions of the 4m condensates obtain 
nonzero values for all five cases listed in Eq. (2), while T x and T Y should belong to the same 
channel. In the case of luld condensates (see Eq. (11 01) ) nonzero values are provided by the 
vector and axial condensates and also by the scalar- vector condensate (T x = 7 , T Y = I) and 

n 

by the axial-tensor condensate (T x = 7 5 7 A , T Y = o^ u ). As noted in [13| . the latter contribution 
was overlooked in [9]. The matrix element Rxy can obtain a nonzero value also in the axial- 
pseudoscalar channel with Tx = 7o7s, Ty =75- However, one can see that the corresponding 
contribution is proportional to q\e XalS0 ■ a a p^. It does not contribute to the structures of the 
polarization operator proportional to q, P or /. 

Note that in the scalar and pseudoscalar cases (r^ = T Y = I, 75) the four-quark condensates 
are determined by one function of density in each of the channels 

H xx = a Xuu (p), R xx = a Xud (p). (11) 
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The same refers to the vector-scalar and axial-tensor condensates (in the latter case Tx = 7a 7s, 
Ty = <Jp,v and a m is proportional to the asymmetric tensor eo^x). The vector condensate with 
Tx = 7 M , T Y = 7^ can be presented in terms of two functions of p 

H X y = a m {p)g^u + b m (p)^f , (12) 

with similar equation for Rxy- The same presentation can be written for the axial condensate. 
Also, there are two functions in the tensor case with T x = o~^ u , IV = cr a p, i-e. 

H X y = a m (p)(g fia g 

) + (13) 



+ b m (py 



m? 



2. Gas approximation 

In the gas approximation the matrix elements in nuclear matter are assumed to be 

(M\X\M) = (0\X\0) + p p (p\X\p) + p n (n\X\n) (14) 

for any operator X. Here \p) and \n) denote the unpolarized nucleons in QCD vacuum, p p and 
p n are the corresponding densities. In the first term on the RHS of Eq. (1141 all four quark 
operators act on the vacuum state. In the gas approximation 

H X y = h p XY p p + h XY p n ; R XY = r XY p p + r XY p n , (15) 

with the matrix elements for unpolarized nucleons 

h N xY = (N\T X y\N); r N XY = ( N\T XY \N) . (16) 

In the gas approximation the functions a m (p) and b m (p) introduced in the previous Subsub- 
section - see Eqs. (fTT|) - (fT5j) are proportional to p. Thus the ratios a m (p)/p and b m (p)/p do 
not depend on p and can be expressed in terms of the nucleon matrix elements hxY an d txy 
introduced in Eq. fflBT) . 

In the gas approximation we can write Eq. (EJ) in the form 

(n m ) 49 - (n ) 49 = ( a4 + A p( ^$- 4 + A 1 ™ 4) ^ p , (17) 

V q m z q z m q z I {27T) Z 
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with 

where the RHS is summed over the flavors /i,/2 and the Lorentz structures X,Y. Here we 
used the standard notation a = — (27r) 2 (0|gg|0) (a > 0). Note that this is just a convenient 
scale for presentation of the results. It does not reflect the chiral properties of the nucleon. The 
coefficients A 1 contain also dependence on the asymmetry parameter 

a = {p n - P P ) J {Pn + p P )- (18) 
Since in the case of 4w quarks only structures with X = Y contribute, we put A^ x = A™. 

C. Three types of terms 

Now we turn to the models in which the nucleon is considered as a composition of the 
constituent quark and the pion cloud. The nucleon vector of state is thus \N) = \(fio,ir), with 
\<po) the state of three valence quarks. Note that all ingredients are in the QCD vacuum. If the 
pion cloud is treated perturbatively, the nucleon expectation value of the operator defined by 
Eq. (j3J) can be written as {9} 

(N\T x f\N) = (MT&lto) + (0o|^/|0o,vr)(vr|T^ 2 |7r)(0 o ,7r|i7 / |0o) + (N\T x f\N) INTERF . 

(19) 

Here Hi is the interaction between the constituent quarks and the pions. In the first term on 
the RHS of Eq. ( Tl9|) all four operators act on the constituent quarks. In the second term the 
four quark operators act on pions - Figs. 1 and 2. 

In the last term two quark operators act on the constituent quarks. Two other ones act on 
the pions. This contact interference is illustrated by Fig. 3a. In the vertex interference shown in 
Fig. 3b the 7rQ interaction Hi adds a pion to the constituent quark. In the four-quark operator 
two operators act on the constituent quarks while the two other ones annihilate the pion. In 
the case of 4m operators we can write 

(N\TT\N) INTERF = ^ ((0 o |^/|0o,7r)(vr|txk)(0o|tx|0o)(0o,7T|^/|0o) + (20) 
+ (Mtx\<M(0\tx\n)(<fo,MHi\<fo) + (0o|^/|0o,vr)(7r|t x |O)(0 o |tx|0o)), 
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with tx = uTxI c u. Since the vacuum and the pion states are colorless, only the first term on 
the RHS of Eq. (3) contributes. The first term on the RHS has a nonzero value in the scalar 
case = /. The two next terms obtain nonzero values in the axial and pseudoscalar cases 
with the matrix elements (Oltxl 71 ") determined by PCAC. Similar equation can be written for 
2u2d condensate. 



D. Contributions to the polarization operator 

Here we consider only those of the four-quark condensates which contribute to the polariza- 
tion operator. 



1. Sea quarks 

In the gas approximation the sea quarks are contained in the pion clouds of separate nucleons. 
If the pions are treated perturbatively, the four quark operators should act on the same pion. 

This requires calculation of the matrix elements (7T a \T X Y \7i a ) with T x p defined by Eqs. 
(pQ) and ([3]), here a denotes the pion isotopic states. These expectation values can be expressed 
in terms of the vacuum expectation values (0|Tj^ 2 |0) by means of the current algebra tech- 
nique 20]. If the factorization hypothesis for the latter expectation values is assumed [jjj], one 



finds that they obtain nonzero values only for X = Y and 

J2{^(^\TTxK)+r x (^\T x d x \^)) = (21) 

x 

for any isotope index a. 

Thus using Eq. ( l2Tj) we find that the contribution of the sea quarks vanishes. 



2. Contributions of the valence quarks. General equations 

The valence quarks are assumed to be represented by the constituent quarks, i.e. by the 
massive particles moving in a certain self-consistent effective field. The constituent quarks are 
assumed to be described by the single-particle functions, which are the solutions of the wave 
equation. In the general case the constituent quarks are relativistic particles, described by the 
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wave function 

^(r) = u q f q (r) (22) 
with / a certain function of r = |r| while in 

u q = f X " ) (23) 

Xg is the two-component spinor with (3 q a certain function of r. In the nonrelativistic models 

p q = o. 

The QCD operator q(r) acts on the valence quark \Q) placed into the QCD vacuum: 

q a \Q a (r)) = Mr)\Q a )+q%r)\0), (24) 

with ip q {r) defined by Eq. ( 1221) . The product of the wave functions describes thus the excess of 
the quark density inside the consistent quark over that of vacuum 

Mr)Mr) = (Q\q(r)q(r)\Q)-(0\qq\0) (25) 

(recall that qq = Sg a g a ). 

3. Nucleon matrix elements 
Here we present the model-independent relations for the nucleon matrix elements. The 



nucleon matrix element for any quark operator X(r) can be written as 21] 



(N\X\N) = (N\ J d 3 r(X(r) - (0|X|0)J|JV> , (26) 

thus expressing the excess of the quark density inside the nucleon over that in the QCD vacuum. 
For the scalar two-quark operator this leads to 



(N\qq\N) = (N\ J d s r[q(r)q(r) - (0\qq\0) J\N) , 
while for the other operators qTxq we have just 

(iV|gr x g|iV) = (N\ J d 3 rq(r)T x q(r)\N) , 
with most of these expectation values vanishing for the unpolarized nucleons. 
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For the four-quark operators qT x qq^Y<l with T x = I , T Y ^ I, i.e. for the operators 
X = qqqTyq one can write Eq. ( 1261) in another form 

(N\qqqT Y q\N) = 2 (0\qq\0) (N\qT Y q\N) + 

+ (N\ J d 3 r(q{r)q{r) - (0|gg|0))g(r)r y q(r)\N) , (27) 

while for the four-quark scalar condensate 

{N\qqqq\N} = 2 (0| qq |0) (N\ qq \N) + 

+ (N\ I d 3 r(q{r) q{r) - (0| qq |0)) V) • (28) 



For the other four-quark condensates qTxqqTyq with fx(Y) = Tx(y) or fx(V") = Tx(y)A a (see 
Eq. (j3J)) we have just 

(AT|gf x ggf y g|A^) = (JV| J Srq(r)t x q(r)q(r)t Y q(r)\N) . (29) 

We shall see that Eqs. ( 1271) . ( 1281) will be useful in the calculation of the contribution of valence 
quarks. 



4- Factorized and internal terms 

In our QCD sum rules analysis there are the factorized terms in the cases of the scalar-scalar 
4u condensate and of the vector-scalar 2u2d condensate. Two quark operators act on the QCD 
vacuum while two other ones act on the valence quarks of the nucleon. We can write 

(N\uI c uuI c u\N) va i = (N\uI c uuI c u\N) f act + (N\uI c uuI c u\N) int , 
(N\u l0 I c udI c d\N) val = (N\u l0 I c udI c d\N) fact + (N\uj I c udI c d\N) int (30) 

- see Fig. 4. Here the lower index vol denotes that we consider the contribution of the valence 
quarks. Index int shows that all four quark operators act on the valence quarks. The scalar 
factorized term is 

(N\uI c uuI c u\N) fact = 2(0|iZ«|0)J, (31) 

where 

3 = 1 d 3 r^ u (r)^ u (r), (32) 
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is just the contribution of one valence u quark to the expectation value (p\uu\p). Another 
factorized term is 

(N\u l0 I c udd\N) fact = (0\dd\0) (N\uj I c u\N) = (0\dd\0) ■ n u . (33) 

Here n u denotes the number of u quarks in the nucleon. Note that the value of the condensate 
(N\ify I c udd\N) f act does not depend on the particular model of the nucleon. 
For the other condensates 

(N\ut x uut Y u\N) val = (N\uf x uuf Y u\N) int . (34) 

The internal terms vanish if a 4m operator is averaged over neutron 9|. For the proton 

(p\T x u x \p) int = 2 J d 3 rt uu r«rf V uu , (35) 

with the wave function of the two constituent u quarks in the proton, the upper indices 
(1), (2) label the two u quarks. Since the asymmetry of the function is provided by that 
of the colors, the spin function is symmetric and the total spin of the two u quarks is 5* = 1. 
Thus we must put (cr^ • cr$) = 1 in the averaging over the spin states in Eq. ( |35i) . For the 
2u2d condensates 

(iV|T! d y |AO ini = Yl i d^ ud T^T^ ud , (36) 

q J 

with the two-particle wave function of a system of u and d constituent quarks. The sum 
is taken over all ud states in the nucleon. Since the spin of the nucleon is 1/2 while the sum of 
the spins of the two u quarks in the proton (and of the two d quarks in the neutron) is 1, we 
must put ^2J(t u ■ o~d) = ~ 4 in the averaging over the spin states in Eq. (|36|) . 



5. Contact and vertex interference 

Start with the 4m condensate. The first term on the RHS of Eq. (1201) describes the contact 
interference, in which two operators act on the constituent quarks while the two other ones act 
inside the neutral pions. There is a nonzero expectation value only for the scalar condensate 

2ml 



(ir\qq\n) 
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with m^, m u and the masses of pion and of the light quarks. 

Similar situation takes place for the luld condensate. Now the charged pions also con- 
tribute. Some of the vector condensates (e.g. (7r°|<i7ott|7r)) have nonzero values. However, such 
expectation values do not contain large factor m n / {m u + mj. Thus their contribution is much 
smaller than that of the scalar term. 

In the second term on the RHS of Eq. ( 1201 one of the vertices of interaction between the 
consistent U quark and neutral pion is replaced by the four-quark condensate. This will be 
referred to as the vertex interference. Nonzero contributions are provided by the axial and 
pseudoscalar condensates. 

In the case of the 2u2d condensates Eq. (7) presents the polarization operator in terms of the 
four-quark operators uY xudTyd. However in the calculations of vertex interference involving 
the charged pions 7r =1= it is more convenient to present the polarization operators in terms of 
the operators uT xddTyu. This is because in the axial and pseudoscalar channels the operators 
uT x d and dT Y u (T X (y) = 7^75, 7s) connect pions and vacuum by the PCAC relations. The 
term J2xy tXY Rxy/<1 2 on the RHS of Eq. (7) is replaced by (see Appendix A) 

^ X,Y 

with 

R XY = (M\f x d Y \M), (38) 
where T ud is written in terms of operators uT xddTyu, and 

i XY = ^757^*7 W y 7 ,75 ■ (39) 

III. CONTRIBUTIONS OF THE FOUR-QUARK CONDENSATES IN PCQM 

A. Wave functions 

In the PCQM the single-particle wave functions defined by Eq. (1221) are 

fu(r) = f d (r) = f(r) = iVe"^ 2 * 2 (40) 

with determined by the normalization condition jip q (r) , yoijj q (r)d 3 r = 1. The two particle 
wave functions take the form of the products of the wave functions given by Eqs. ( 1221) and ( 1401 
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with the spinors xi,2 composing the corresponding spin state. The functions /3 q are replaced by 
constants 

&u = Pd = P = 0.39, (41) 

fitted to reproduce the value of the axial coupling constant. The parameter R is related to the 
contribution of the valence quark core (r 2 ore ) to the proton charge radius (r 2 ) [11 ]. 



In the PCQM calculations of the proton charge radius [22| the value 

R = (0.618 ± 0.006) fm 



provides the exact fit for the nowadays value (r 2 ) 1 / 2 = (0.8768 ± 0.0069)fm 23] 

A characteristic size of the internal expectation values is 

1 1 1 



f\r)d A r 



(2tt) 3 / 2 ' R 3 (1 + 3/2/3 2 ) 2 



with the numerical value for R = 0.618 fm 

Af 2 = 1.37- 10~ 3 GeV 3 

Here 



0.098 e 3 . 



4 = -(0\qq\0), = 241 MeV, 

thus the second factor on the RHS of Eq.(17) is a/(2ii) 2 = Eq. 



(42) 



(43) 



(44) 



(45) 



B. Contributions of the valence quarks 



1. Factorized terms 



Employing the PCQM value J = 0.54 we find for the factorized term provided by Eq. ( 1311 ) 

A% = 1.08 (l - | ) , (46) 

with the asymmetry parameter a defined by Eq. fflBl) . For another factorized term expressed 
by Eq. (|33|) we obtain 

A^ s = 2.00 (l - |) ■ (47) 

Inclusion of nonlocality of the vector condensate and of some other small corrections [9| changes 
the last equation to (see Appendix B) 



(4* 
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2. Internal terms 

The internal contributions of the valence quarks to the parameters A^' 1 defined by Eq. (p3 
can be obtained by straightforward calculations (see Appendix C). We find the numerical values 
for the 4u case 9] 

A% = -0.05 + 0.04a; A q Ps = -0.01 + 0.01a; Ay = 0.01(1 - a); 

A\ = -0.03(1 - a); A% = 0.02(1 - a) , (49) 

and 

A p = A p Ps = ty Ay = -0.13(1 - a); A P A = 0.03(1 - a); A£ = -0.09(1 - a), (50) 

while A 7 = for all these structures. 

In the case of the 2u2d condensate the nonvanishing contributions are 

A v = -0.80; A\ = -0.49, (51) 

while 

A v = -0.41 , A P A = -0.25, (52) 

and 

A vs = -0.34, A AT = 0.77. (53) 
Thus the total contribution of the internal terms is 

A q = -1.35 + 0.05a, A p = -0.85 + 0.19a, (54) 

A 1 = 0.43. 



3. Total contribution of the valence quarks 

Summing up the contributions provided by Eqs. (|46j) . (148]) and (|54j) we find for the total 
contribution of the valence quarks 

A q = -0.27 - 0.31a, A p = -0.85 + 0.19a, (55) 

A 7 = 1.57 - 0.38a. 

One can see that there is a strong cancelation between the factorized and internal terms in the 
structure A q . Internal terms provide a small correction to the structure A 1 . 
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C. Interference terms 

We now turn to the case in which one of the operators qT x,Y q acts on the constituent quark 
while the other one acts on the pion. In the PCQM the pions manifest themselves in the ir a Q 
self-energy loops E a and in the corresponding exchange terms I a - see Fig. 3. They are related 



as |a] 



I a = (56) 



with the same relation for the contributions of the four-quark condensates. 
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For the contact interference described by the first term on the RHS of Eq. (120]) we find [9] 
(see Appendix D) 

A q = -0.02 + 0.01a; A 1 — —0.13 + 0.06a. (57) 



For the vertex interference, described by the second and third terms on the RHS of Eq. (1201) 
we obtain in the case of 4m condensate 

(Q\TTx\Q,7r) = ^{Q\uT x u\Q){U\uT x u\ir). (58) 

The matrix elements {0\qT x q\iT) obtain nonzero values for the axial and pseudoscalar op- 
erators. Similar equation can be written for the 2u2d condensate. Since the ttQQ vertex is 
proportional to the product (er • k), where k is the pion momentum, the vertex interference can 
take place if the product {Q\qY x q\n) {Q\qYyq\Q) also has this structure. 

Using the definitions 

fn = 97^75? ; ji = u-i^d ; j~ = d^ 5 u , (59) 
one can write the PCAC relations 

(0b>°(A:)) = <M.. (0\j d y(k)} = -iUk^ (0|j;|vr + (fc)) = iyfif^, (60) 

and 

(0|u7b«|7r°(*)> = (0\d l5 d\7i°(k)) = i f -^-- (0|rf 7 5«|vr + (A:)) = i f ^ , (61) 

2m u 2m d m u + m d 

which will be used in calculations of the four-quark condensates. 
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If both Tx and Ty are the axial matrices, the contribution of the four-quark condensate can 
be expressed in terms of the self-energy loops of the constituent quarks (see Appendix D). We 
find 

A q = -0.65 + 0.02a; A p = -0.04 - 0.02a. (62) 

The large value of A q is provided mostly by the charged pions. 
If T x = T Y = 7 5 

A q = 0.84 -0.07a, (63) 

also with the main contribution provided by the charged pions. 

There is also a contribution of the axial T x and tensor T Y . It leads to 

A 1 = 0.36 + 0.04a. (64) 

Thus for the sum of the interference terms - see Eqs. (57), ( |62l) -(l64l ) ) 

A q = 0.17 - 0.04a , A p = -0.03 + 0.02a , (65) 

A 1 = 0.23 + 0.10a. 

D. Total contribution 

The sum of the values given by Eqs. (55) and ( |65l) provides for the coefficients on the RHS 
of Eq. (HZD 

A q = -0.10 - 0.34a, A p = -0.88 + 0.17a, (66) 
A 1 = 1.80 - 0.28a. 

In Table 1, we present the contributions of various Lorentz structures to the matrix elements 
Hxx defined by Eq. (E]) and to the parameters A 1 (Eq. ( flTl) ) for the Au condensates. In Table 2, 
we show similar data for the 2u2d condensates. In the latter case, all the contributions except 
the vertex interference terms involving the charged pions are calculated in terms of operators 
uY xudTyd. The interference terms involving the charged pions are calculated in terms of the 
operators uT xddTyu. Note that the data in the Tables 1,2 contain the contributions of the 
sea quarks for each Lorentz structure. They cancel in the sum over all Lorentz structures - see 
Eq. (21). 
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E. Calculations in the nonrelativistic quark model 

It is instructive to compare the results with those in the old nonrelativistic quark model 
(NRQM). The nucleon is treated as a system of three nonrelativistic constituent quarks in a 
self-consistent central field. The wave functions can be taken in the form provided by Eqs. 
( l22l) . ( )23l) and (HDjl with f3 q = 0. We must equal (r 2 ore ) to the proton charge radius. This leads 
to R = 0.72 fm. The NRQM values for the symmetric matter are 

A q = 1.15, A p = -0.54, A 1 = 1.45. (67) 

The strong deviation between the PCQM and NRQM values for A q is due to the difference 
between the values of the contribution of the valence quarks to the condensate (N\uu\N) - see 
Eq. (132 p . Due to relativistic reduction the PCQM value is about twice smaller than the NRQM 
one J — 1. 



IV. SOLUTIONS OF THE SUM RULES EQUATIONS 



In this Section, we present equations for symmetric nuclear matter (a = 0). Following |4j. Il7] 
we consider the dispersion relations Eq. (4) at fixed value of the relativistic pair energy 

s = (P + q) 2 = 4m 2 . 

This enables to separate the singularities connected with the nucleon in the matter from those 
of the matter itself. Now we clarify the terms on the LHS of Eqs. (jHJ). The contributions of 
continuum are expressed in terms of the functions E{(x) with x = W 2 /M 2 

E (x) = 1 - e~ x - E x {x) = 1 - (1 + x)e~ x ; E 2 (x) = 1 - (1 + x + x 2 /2)e~ x . (68) 

The radiative corrections of the order a s and a s In M 2 / A 2 to the condensates of dimension n con- 



tained in the structures A n , B n and P n - Eq. ([6]) [24| are expressed by the factors fA(B,p)n with 
?A(B,p)n — 1 if the corrections are not included. We assume the numerical value A = 230 MeV. 
In the one-loop approximation this corresponds to a s (lGeV 2 ) = 0.475. The contribution of 
the lowest dimension is [15[ 

Aq = M 6 E 2 (W 2 /M 2 )f A0 , (69) 



19 



corresponding to the loop of three free quarks with inclusion of the radiative corrections. 
Next come the terms, containing the in-medium condensates of the lowest dimension d = 3. 



A 3 = _ 8l As-^)E a (wyAfi)-M^(WV^) 

3m 

B 3 = -A7i 2 M 4 E 1 (W 2 /M 2 )K(p)f B3 ; 

3?TT 2 

P 3 = — M 4 E 1 (^ 2 /M 2 )t;(p)fp3, (70) 



3 

with the vector and scalar condensates 

v(p) = (M\Y^fa tf\M) = v N p; v N = 3, (71) 

i 

and 

k( P ) = (M\Y^MM), (72) 

i 

with the sum over the light quarks. 

The nonlocality of the vector condensate and the gluon condensate compose the contributions 
of dimension d = 4. They are numerically small, thus we do not display them here. However, 
we take them into account in our calculations. The contributions of dimension d = 6 come 
mostly from the four-quark condensates. In vacuum (p = 0), the only contribution to the 
polarization operator is (see Eq. <^§) 

(n ) 4g = -- ^ ' 2 ' ' qr Ae , (73) 
o q 



providing 



4 

A 6 = -a 2 r A6 . (74) 



Recall that a = — (27r) 2 (0|gg|0) > 0. In nuclear matter we must add the contributions of the 
nucleons. 

A. Gas approximation 

The contribution of the density- dependent part of the four-quark condensates is given by 
Eq. ( fl4l) . Their contributions to the LHS of Eqs. (JSJ are 

A% = -8n 2 aA g p ; B 9 6 = -12n 2 am A 1 p; P% = -12n 2 am A p p , (75) 
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with the parameters A 1 given by Eq. (166 p . The upper index g denotes that the values are 
calculated in the gas approximation. 

In the gas approximation the scalar condensate determined by Eq. f!72|) can be written as 



«(P) = (0| ^2^\0)+k p ) k p = k nP ; k n = (N\J2<fl l \N). (76) 



Parameter can be expressed in terms of the pion-nucleon sigma term 25 ] 



k n = ; , (77) 



while a n N is related to observables 
term are somewhat controversial 



26 
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271 ] . The experimental data on the value of the sigma 



291 ]. providing rather an interval of the possible values 

7 < k n < 12. (78) 

We present the density dependence of the nucleon parameters in Fig. 5. The numerical 
results at the phenomenological value of the saturation density po — 0.17 fm -3 are given in 
Table 3. One can see that at = 8 inclusion of the four-quark condensates subtracts about 
100 MeV from the value of the vector self-energy and about 170 MeV from the value of the 
effective nucleon mass. Also the residue of the nucleon pole becomes about 2/3 of its vacuum 
value. The influence of the four-quark condensates on the nucleon parameters is even larger for 

K N = 11. 

Employing Eq. (5) we can write 

VK ! Ci{p-M 2 ,Wl) yi h y 1 Ci{p-M 2 ,Wl) KHJ v ; 

The values which approximate the functions Fyj(M 2 ) can be treated as the solutions of the 
sum rules equations. Replacing the effective threshold W^(p) by its vacuum value Wq we find 
the approximate solution of the sum rules equations 

C p (p;M 2 ,W 2 ) C^M^Wj) 

C»{p;M*,W2) =- S ^ ; £*(p;MW ') =m(p) ' (80) 
which solve the problem of expressing the parameters of the nucleon in nuclear matter in terms 
of the density-dependent QCD condensates. Putting = Wq we find Sy = 253 MeV and 
m* = 599 MeV for p = p and = 8, while the solutions of Eq. ( I79|) are Sy = 217 MeV and 
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m* = 562 MeV - see Table 3. Hence the accuracy of the approximation = Wq for the 
self-energies £y and m*—m is about 15 percent. It becomes less accurate for larger values of 
kn providing the errors of about 20 percent at — H- 

Note that inclusion of the four-quark condensates lead to m* + < m i i- e - the nuclear 
matter is bound. 

The large impact of the four-quark condensates on nucleon parameters is to some extent 
caused by employing the gas approximation. 

B. Beyond the gas approximation 

Now we take into account the contributions to the four-quark condensates, in which two pairs 
of quark operators act on two different nucleons of the matter. We include only the contributions 
containing the large parameter k^. These are the scalar and scalar- vector contributions. 

The contribution proportional to be obtained by replacing the factor ((0\uu\0)) 2 on 

the RHS of Eq. ( 1731) by (k p /2) 2 - Eq. ( 176|) . It comes from the 4u condensate and contributes 
to the q structure. In order to obtain the contribution proportional to the product knVn it is 
sufficient to replace the vacuum expectation value (0|<i<i|0) in Eq. ( 1331) by tt p /2. This contributes 
to the scalar structure / comes from 2u2d condensates. 

The corresponding contributions to the matrix elements defined by Eq. flH]) are 

H S s = ^%P 2 ; Rvs = ^v N n N p 2 . (81) 
The contributions to the LHS of Eq. ([5]) are now 

167T 4 

A 6 = A% + ^-4P 2 ! B 6 = B 9 6 + 8tt a v n k nP 2 ; (82) 

h = Pi- 

Thus the large contribution of the four-quark condensate in the scalar channel A 1 (see 
Eq. ( 16"B"|) ). obtained in the gas approximation becomes about 30 percent smaller due to inclusion 
of the 3-body interactions. 

The density dependence of the nucleon parameters is shown in Figs. 5 and 6. The numerical 
results for p = po are shown in Table 3. One can see that the influence of the four-quark 
condensates on the values of the nucleon effective mass and of the residue becomes weaker, 
than it was in the gas approximation. 
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C. Results in the nonrelativistic quark model 

Now we present the results with contributions of the four-quark condensates calculated in 
framework of the NRQM . We find that for = 8 they differ strongly from those obtained in 
the PCQM - see Table 4. The vector self-energy £y is about 300 MeV, and the effective mass 
m* differs from the vacuum value only by about 100-130 MeV. The nuclear matter appears 
to be unbound since m* + £y > m. This is due to large factorized term in the q structure 
which is not compensated by the internal terms in the case of NRQM. For kn — 11 the scalar 
term increases the value of m — m* and the matter becomes bound. However the values of the 
parameters differ much from those predicted by the PCQM. 



V. SUMMARY 



We calculated the four-quark condensates in nuclear matter in the gas approximation and 
included the most important terms beyond it. This corresponds to inclusion of the two-body 
forces and of the most important three-body forces. We applied the results for solving the QCD 
sum rules equations (Eq. ([5])). 

The calculation of the four-quark condensates requires certain model assumptions on the 
structure of the nucleons which compose the matter. In the broad class of models the nucleons 
are treated as the systems of valence quarks, described by the constituent quarks and the sea 
quarks, contained in pions. We showed that in the models which treat the pions perturbatively, 
there is a remarkable cancelation of the contributions to the sum rules which contain only the 
sea quarks. Among the various contributions which involve the valence quarks one can separate 
the factorized terms, where two operators act on the nucleons of the matter, while two other 
ones act on the QCD vacuum and internal terms in which all four operators act on the valence 
quarks. In the interference terms two of the operators act on the valence quarks. The other 
two ones act on the sea quarks (i.e. on the pions) or connect the pions with vacuum by the 
PCAC relations. 

Calculations, carried out in framework of PCQM showed that the contributions of the four- 
quark condensates to the scalar channel A 1 can be viewed as dominated by the factorized terms. 
The contribution to the vector structure A p is almost totaly due to the internal terms. The 
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vector structure A q is a result of strong compensation of the factorized term given by Eq. (1461) 
by the internal and interference terms. 

Solving the sum rules we find the nucleon self-energies presented by Eq. (179 p . Approximate 
solution presented by Eq. (|80|) solves the problem of expressing the nucleon self-energies in 
terms of the QCD condensates. 

The values of the nucleon parameters at the phenomenological value of the saturation density 
are presented in Table 3. One can see that the large reduction of the nucleon effective mass 
m* and of the nucleon residue is to some extent an artifact of the gas approximation. The 
effect becomes less pronounced when the most important terms beyond the gas approximation 
are taken into account. 

Anyway, the nuclear matter at the phenomenological saturation value appeared to be bound. 
The values of the four-quark condensates are consistent with the hypothesis of the convergence 
of the power series in 1 /q 2 of the polarization operator. Uncertainties in the experimental data 
on the nucleon matrix element k n = (N\uu + dd\N) (Eq. (I7H!) ) lead to uncertainties in the 
quantitative results. 

The density dependence of nucleon parameters found by employing the PCQM are shown 
in Figs. 5,6. Also, in Fig. 7, we compare the results for two values of kn- 

It is instructive to have a look at the results obtained in the nonrelativistic quark model. 
The contribution of the four-quark condensates to the coefficient A q differers strongly from that 
in PCQM - Eq. (16"T|) since there is no relativistic reduction of the contribution of the valence 
quarks to the expectation value (N\uu\N) this time. The results for the nucleon parameters 
presented in Table 4 contradict to the nuclear phenomenology (the matter appears to be not 
bound) for k>n = 8. It is bound for = 11 but the quantitative results differ from those, 
obtained in the PCQM. 

The results of further investigation with more complete inclusion of the many-body forces 
will be published elsewhere. The authors acknowledge the partial support by the RFBR grants 
11-02-00120 and 12-02-00158 and by the grant RSGSS -65751.2010.2. 
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Appendix A. 

Here we calculate the contribution of the four-quark condensates to the polarization operator 

n m (g 2 ) = i [ d^xe^(M\Tj(x)](0)\M). (A.l) 



We employ the current suggested in [18]. It is 



h (x) = e abc (<C 7 ^)(757^ c ) 7 , (A.2) 

where u and d are the quark operators, a, b, c are the color indices, C is the charge-conjugation 
matrix (C T = C^ 1 = — C) and the upper index T denotes the transposition. Since we want to 
calculate U m (q 2 ) as a power series of q~ 2 , we need expansion of the function 

jr c ( x ) = Tj 7 (x)j c (0) = -^ abc £ de ^^(x)(C 7 , t )^4(x)(757 /i ) 7 «5^(x) x 

x dl(o)(Y^U%(oWC) vP mt ■ (A.3) 

in powers of x 2 . For a product of the two quark operators we can write 

95U(*)?U°) = G Jw^-^E^°) r 4(°) r J^ - 

X 



with the fundamental matrices defined by Eq. (2), f± are the quark flavors, while G ab (x) 
G(x)5 ab where 



G(x) = (A.5) 



is the propagator of a free massless fermion. 

In the case of 4m condensates, we describe the product d c s (x)dl(0) in Eq. (1A.3j) by the first 
term on the RHS of flA.4j) . while the other products of the field operators are given by the two 



next terms on the RHS of that equation. Since 



d 4 xG(i)e i(9l) = (A.6) 
q 2 



we come to Eq. (9) after doing some matrix algebra. We employed that CT^C = ±T X with 
the plus sign in the vector and tensor cases and the minus sign in other cases. 
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In the case of 2u2d condensates the product of two u quark operators (e.g. u^(x)u^(0)) is 
described by the propagator G ab (x) given by Eq. ( 1A.5I) . Employing Eq. ( 1A.4[) we describe the 
product of the operators uudd in terms of the operators uTxU ■ dTyd, coming to Eq. (10). An 
alternative presentation in terms of the operators uZxd ■ dZyu can be obtained by employing 
the commutation relations for the quark operators and presenting the products ud and du 
employing Eq. ( 1A.4I) . This leads to Eqs. ( l3Tj) -( l39l) . The two presentations are tied by the Fierz 
transform. 



Appendix B. 

The general expression for the contribution of the factorized vector-scalar term (Eq.(27)) to 
the LHS of Eq. (4) is [l9| 

rw,) = -H<?!fW/^(«. SW ,«*-) ft , (B.i) 

with 6{x) = 6 u (x) + 6 d (x), while 

8l( x ) = (p\<l(Qhn<l( x )\p) = —Vai^iX^ + ix^mipbdpx)^ 2 ), (B.2) 

where 

q(x) = q(0) + x a D a q(0) + ... (B.3) 
In the leading order we can put x 2 = in Eq. flB.2p and define HQ 



Va{b ){{px)) = ip a{b) {(px,0)) = / daf a(b) (a)e-^ a . (B.4) 



with f a {<y) = f{ct) the proton deep inelastic structure functions 30j. The moments of the 
functions fb{oe) can be expressed in the terms of the moments of the functions f{ot) 171 ]. 
Including the two lowest moments (a"" 1 /) we find 

2(0|dd|0) (2{pq) lex A(pq) 
Tq 

Using the normalization condition (/) = 3 and the numerical value (otf) = 0.45 [3jJ we come 
to Ay S = 1.30 (1 — a/3) . Inclusion of the lowest correction to the vacuum expectation value 
(0|dd|0) [19( leads to Eq. flS}. 



2(0\dd\0) f2{pq) /n A{pq) \ 

n(g)/act = =-3 — (/> (af) -m(af))p p . (B.5) 

3a V m m / 
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Appendix C. 

Here we give several examples of calculation of the internal terms - Eqs. (1351) and fl36|) . For 
the vector and axial terms the matrix elements for the 4m condensate defined by Eq. (fl6|) is 

h% = + (C.l) 

with fi(v) the indices of the vector and axial matrices 7 M and 7^75. There is similar equation 
for the 2u2d condensate described by the matrix element r^ u . 



In the vector case one can calculate for the time components 



Ko = a v + b v = I d 3 xf(x) ( 1 - /3 2 -^) , (C.2) 



2 x 2 



l hUij = -a v = ^ f d 3 xf(x)^w, (C.3) 



and for the space components 

3 13 3 3 J J v ' R 2 

with 

w = 2(x\ct^. ( t^\x}. (C.4) 

Here x is the spin wave function of the nucleon, while 1 and 2 label the two u quarks. Since 
the total spin of two u quarks is S = 1, we find w = 2. 

Similar equations can be written for the 2u2d condensate. However this time 



w 



J2(x\<ru-<r d \x), (C5) 



with the sum over all quarks of the nucleon. For the proton w = (x|(<Xu + cru"') ■ o"d\x)- Since 
the spin of the proton is S = 1/2 we can write l/4(crl 1} + crl 2) + cr d ) 2 = 3/4. Employing 
((Tu^ ■ CTu) = 1 we find w = —4. In the same way we find w = —4 for the neutron as well. 
This provides for the 4m condensate 



a v = _ V 2 /3 2 = -0.20^; b v = 2Af 2 (l + —f3 2 + -/3 4 ) = 0.26 eg , (C.6) 

3 V 6 16 / 

while for 2u2d condensate 

a v = V 2 /? 2 = 0.044 b v = 2J\f 2 (l + -/3 2 + -/3 4 ) = 0.21 eg . (C.7) 

3 V 6 16 / 
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n 

One can see that in axial case the time components turn to zero [9j. Thus 

a A + b A = 0, (C.8) 

Employing the values of fi XY and t xy (see Eqs. (9) and (10)) we find the contributions 
presented by Eqs. (B9l-(]52l). 

In the axial-tensor case T x = 7s7 A (r^ = 7a7s)? T y = ^ v \ we obtain for the ingredients of 
Eq.© 

Tr( 7 a g 7 /3 r x ) = Tr(7 a g 7 /3 7 5 7 A ) = Tr(g 7 a 7 /3 7 A 7 5 ) = Aie Ka ^ x q K , (C.9) 

while 

757«r y 7/375 = iflf W *' ' g W '(V«J^ - 9n'p9u' a ) , (C.10) 

and thus 

t at„u\ = 2q p£ P^ x . (C.ll) 
Employing Eq. (1231) for u q we find that the axial and tensor matrix elements are 

^7075^ = 0; cu q 'ya 5 uj q = a { + (3 2 ( crx )^( <TX ) ? (i = 1,2,3), (C.12) 

Xj- _ / 2 (crx)(T/ c (crx) \ . 
u q a 0j uj q = -2(3—; u q cr i:j uj q = e ijk \cr k - p — J, [% = 1,2,3), (C.13) 

where is the standard three-dimensional asymmetric tensor. 

Thus the indices /i, v and A on the RHS of Eq. flC.llj) should correspond to the space com- 
ponents, while p = 0. Hence we can write 

t at,uX = 2qo£ oi jk = _ 2 q s 0ljk , (C.14) 

with i,j, k - the space components, and £0123 = 1- Hence we can put £ y& = £ij k - 
Thus 

2qo ^>I CT « • "tlx) J d 3 xf 2 (x) (l - ^ . (C.15) 
Hence 

This leads to Eq.(53) 

A AT = 0.77. (C.17) 



ATfivX AT 
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Appendix D. 

Here we calculate some of the interference terms. The contact interference and the vertex 
interference in the axial case are described in in details. Consider here the pseudoscalar 
vertex interference. 

Start with the case of the Au condensate. Only the neutral pions contribute to the interfer- 
ence terms. For the matrix element h = hp s defined by Eq. (1161) we obtain {9] 



8 1 Q 

h = ~ J°n u . (D.l) 

Here J° is the result of replacement of one of the UUn vertices in the self-energy loop of the 
U quark 

S° = ~ J d 3 xd 3 zg*(x)S(x)D(x - z)S{z)g{z) , (D.2) 

by the factor —if^m^/2m u - see Eq. (l6T|) . i.e. 

77?^ 1 f 

jo = _ 2L/ i=- d 3 xd 3 zg*(x)S(x)D(x - z)g(z) . (D.3) 
2m u 2 J 

In these equations 

g(x) = i>(x)i^(x) , 



D is the pion propagator, while 11] 



y/20 f ni-yi + r 

1 dy V - ° ' e-y. (D.4) 



Six) - + ^ 

b[X) ~ 2PR + 2& X 

is the scalar field. 

Using the explicit expressions for the functions ip{x) we obtain 

2 

4tt 2 (1 + 3/2/? 2 ) 2 J R 2 J ^ y 2 T^ 

Direct computation provides J° = 5.20 • 10~ 2 e 3 . 

Thus we obtain h = 0.27n u . Since fi Ps = q/2 (Eq. (9)), we find the contribution 

A q = 0.20 - 0.07a. (D.5) 

In the case of 2u2d condensates we write the contribution in terms of the operators uj^ddj^u 
- see Sec. 2.4.4. Defining the matrix element in Eq. fl38|) as Rp s = (M\Tp d s \M) = rp. We obtain 

r=~J C (n u + n d ); J c = , (D.6) 

3 9 m u + m d 
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with the numerical value J c = 7.56 • 10 -2 £q , leading to 

A Q = 0.64, (D.7) 

with the total contribution of the pseudoscalar terms (Eq. (63)) 

A q = 0.84 -0.07a. (D.8) 

Turn now to the axial-tensor interference. Only the luld operators contribute. Start with 
the case of neutral pions tt°. The ttqDD vertex of the D quark self-energy loop can be replaced 
by the axial-tensor condensate 

tx^u = \($\3 u x \AWa^d\D) ■ (D.9) 

Since only the space components contribute to polarization operator (Eq. (1C.12j) ). we obtain 

2 

kij = ^fn£ijr(D\dk e a T d\D) , (D.10) 

with the three-dimensional indices I corresponding to the four-dimensional ones \[iv. Using 
Eq. ( 1C.12I) and employing Sij r £ijt = 25 r e we obtain 



Trfttij = -^Erijttij = -Uu{D\d(k-(T)d\D). (D.ll) 

ij 

This leads to 

h = ffj AT n d] J AT = ±1 = 2.89 • 10- a eg , (D.12) 

leading to 

A 1 = 0.12 + 0.04a. (D.13) 
Proceeding in the same way in the case of charged pions and using Eqs. (137"]) -( 139]) we find 

A 1 = - — J AT (n u + n d ), 

providing 

A 1 = 0.24. (D.14) 
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TABLE I: Contributions of various Lorentz structures to the matrix elements Hxx defined by Eq. ([8]) 
and to the parameters A 1 (Eq. ()17p ) for the 4u condensates. 



X 


a m /(pe ) 


W(P £ o) 




A p 


A 1 


s 


-6.34 + 0.60a 




3.17- 0.30a 








Ps 


1.74 -0.12a 




-0.87 - 0.06a 








V 


-0.55 + 0.01a 


0.13 - 0.13a 


0.55 - 0.01a 


-0.13 + 0.13a 





A 


0.33 + 0.06a 


0.07 - 0.05a 


0.47 + 0.05a 


0.07 - 0.05a 





T 


0.02 - 0.02a 


-0.04 + 0.04a 


0.02 - 0.02a 


-0.09 + 0.09a 






TABLE II: Contributions of various Lorentz structures to the matrix elements Hxy and Rxy de- 
fined by Eq. ([8]) and to the parameters A 1 (Eq. (]17p ) for the 2u2d condensates. All the contributions 
except the vertex interference terms involving the charged pions are calculated in terms of operators 
uT xudTyd. They are given in the lines above the horizontal line. The vertex interference terms 
involving the charged pions are calculated in terms of the operators uT xddTyu. Their contributions 
are presented in the lines below the horizontal line. 



XY 


a m /{p£o) 


W(/° £ o) 


A q 


A p 


A 1 


vv 


0.58 


0.21 


-6.22 


-0.41 





AA 


-0.22 


-0.16 


2.60 


-0.32 





vs 


-0.33 + 0.16a 










0.67 - 0.32a 


AT 


-0.89 - 0.04a 










0.89 + 0.04a 


PsPs 


0.64 




0.64 








AA 


0.23 





-0.46 








AT 


0.24 










0.24 
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TABLE III: Nucleoli parameters in symmetric nuclear matter at the phenomenological saturation 
value of the nucleon density. The first line shows the vacuum values. The results marked as (*) in the 
second and fifth lines are for the vacuum values of the four-quark condensates. The values (GA) in the 
third and sixth lines correspond to the gas approximation. The values in the fourth and seventh lines 
correspond to inclusion of the three-body interactions, following Eq. (|82j) . The results are presented 
for two values of the parameter kn - Eq. (|78|) 





E v , MeV 


m*, MeV 


X 2 m , GeV 6 


Wl , GeV 2 


p = 





940 


1.94 


2.16 


p = p , K N = 8 (*) 


329 


737 


2.39 


2.33 


P = Po, n N = 8 (GA) 


217 


562 


1.21 


1.62 


P = Po, k n = 8 


206 


602 


1.41 


1.69 


p = po, K N = 11 (*) 


335 


606 


1.91 


2.12 


p = Po, k-n = 11 (GA) 


203 


388 


0.98 


1.52 


P = PO, K N = 11 


180 


445 


1.20 


1.58 



TABLE IV: Nucleon parameters in symmetric nuclear matter at the phenomenological saturation value 
of nucleon density obtained by employing of the nonrelativistic quark model. The values (GA) in the 
first and third lines correspond to the gas approximation. The values in the second and fourth lines 
correspond to inclusion of the three-body interactions, following Eq. (|82p . The results are presented 
for two values of the parameter kn - Eq. (JTSj) 





£„, MeV 


m*, MeV 


A^GeV 6 


^,GeV 2 


p = po, kn = 8 (GA) 


305 


685 


1.60 


1.89 


p = P0, KN = 8 


284 


702 


1.84 


2.00 


p = po, K N = 11 (GA) 


291 


504 


1.11 


1.59 


P = P0, K N = 11 


253 


530 


1.39 


1.71 
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Figure Captions 

FIG. 1. Contribution of the four-quark condensates (dark blob) to the left-hand side of the 
sum rules caused by the valence quarks (dotted lines). The helix line stands for the nucleon 
current. Solid lines denote the quarks of the current. 

FIG. 2. Contribution of the four-quark condensates (dark blob) to the left-hand side of the 
sum rules caused by the sea quarks. Wavy line stands for the pions. The other notations are 
the same as in Fig. 1. The irQQ interaction vertices 1 and 2 may correspond to the same or 
different nucleon quarks. 

FIG. 3. Contribution of the four-quark condensates (dark blob) to the left-hand side of the 
sum rules caused by the interference terms with Fig. 3a and Fig. 3b illustrating the contact and 
vertex interference correspondingly. The notations are the same as in Fig. 2. The vertices 1, 2 
and 3 may belong to the same or different nucleon quarks. 

FIG. 4. Contribution of the four-quark condensates to the left-hand side of the sum rules 
caused by the factorized terms. The small circles denote the QCD vacuum. The dark blob 
stands for the two-quark condensate. 

FIG. 5. Density dependence of the nucleon vector self-energy Sy and of the nucleon effective 
mass m* in the symmetric nuclear matter for kn = 8 corresponding to Eq. (79). The horizontal 
axis corresponds to the density p related to its saturation value po- The dashed lines show the 
results in the gas approximation. The solid lines correspond to inclusion of the most important 
3-body forces. 

FIG. 6. Density dependence of the nucleon residue \ 2 m and of the effective threshold 
for Kn = 8 in the symmetric nuclear matter corresponding to Eq. (79). The horizontal axis 
corresponds to the density p related to its saturation value p . The other notations are the 
same as in Fig. 5. 

FIG. 7. Density dependence of the nucleon vector self-energy Sy and of the nucleon effective 
mass m* in the symmetric nuclear matter corresponding to Eq. (79) with inclusion of the 3-body 
forces. The horizontal axis corresponds to the density p related to its saturation value p . The 
solid and dot-dashed lines show the results for kn = 8 and kn — 11 correspondingly. 
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